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On the central black hole mass in Mkn 501 
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Abstract. We analyse the apparent disagreement between the mass estimates of the central black hole(s) in 
Mkn 501 based on (i) the observations of the host galaxy, (ii) the high energy (HE) emission mechanism, and 
(iii) the modulation of the beamed radiation by a black hole (BH) binary system. While method (i) seems to 
imply a central mass ~ 5 x 1O 8 M0, method (ii) suggests a BH mass less than ~ 6 x 1O 7 M0. We critically 
discuss the estimates inferred from (i) showing that current uncertainties may permit a central mass as low as 
~ (2 — 3) x 10 8 Mq. We demonstrate that in this case the estimates (i) and (ii) might be brought into agreement 
by assuming a binary BH system where the jet dominating the HE emission originates from the less massive 
(secondary) BH as suggested by method (iii). On the other hand, if Mkn 501 has in fact a high central BH mass 
of order 10 9 Mq , a change of fundamental assumptions seems to be required in the context of several HE emission 
models. We show, that in this case a binary scenario following (iii) may be still possible if the jet which dominates 
the emission emerges from the more massive (primary) BH and if the binary evolution passes through phases of 
super- Eddington accretion and/or decreased conversion efficiency. 
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particularly relevant for our understanding of BL Lac ob- 
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. According to the commonly accepted paradigm, Active 
Galactic Nuclei (AGN) are thought to harbour supermas- 
sive black holes (BH) surrounded by geometrically thin 
accretion disks, the latter triggering the formation of rel- 
ativistic jets. Increasing evidence now indicates, however, 
that this picture probably has to be expanded: First, hi- 
erarchical galaxy evolution schemes suggest, that as a 

: result of mergers between galaxies, binary BH systems 
(BBHSs) should be generally expected in the center of 
elliptical galaxies, and hence, for example, in the typi- 
cal hosts of BL Lac objects (cf. Begelman et al. 1980; 



jects such as Mkn 501 (e.g. De Paolis et al. |2002|). 
Mkn 501 (z = 0.034) is one of at least four AGN which 
have been reliab ly detected at TeV energies (see Catanese 
& Weekes 1999 for a review). From optical observations 
its host is known to be the elliptical galaxy UGC 10599 
(cf. Stickel et al. |l993|). As a BL Lac object, Mkn 501 



Kauffmann 1996; Richstone et al. 1998; Milosavljevic & 



Merritt |2001| ; Yu |2002[ ) . Secondly, from a phenomenologi- 
cal point of view, a multitude of observational evidence has 
been plausibly related to the presence of BBHSs in AGN, 
including the observed misalignment, precession and wig- 
gling of extragalactic jets (e.g. Conway & Wrobel 1995 ; 
Kaastra & Roos 1992 ), periodic optical outburst eve nts as 
observed in the q uasar OJ 287 (e.g. Sillan paa et al. 1988 ; 
Lehto & Valtonen |l996| ; Villata et al. |l998[ ) and th e helical 
motion of knots (e.g. Abraham & Carrara 1998; Britzen 



belongs to the blazar class of AGN which are thought to 
have relativistic jets oriented at a small angle to the line 
of sight, so that the observed flux is strongly enhanced 
by relativistic beaming effects. Detailed observations since 
1997 have confirmed the picture of powerful activity in 
Mkn 501 and revealed a wide range of variability pat- 
terns during outbursts (e.g. Prothcroe et al. 1998) includ- 
ing flaring episodes of several days and rapid variability 
on subhour timescale (e.g. Sambruna et al 2000] ). The evi- 
dences for a 23-day periodicity in both the TeV and the X- 
ray lightcurves during the 1997 high state (cf. Hayashida 
et al. |1998| ; Kranich et al. |1999| ; Nishikawa et al. |1999[ ) 



may be counted among the most fascinating features (see 
Kranich 2001 and Kranich et al. 2001 for a recent assess- 



et al. POOID . 

One expects that a binary framework for AGN may be 
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mcnt of significance) and possibly indicates the presence 
of a BBHS in Mkn 501 (cf. Rieger & Mannheim p000| ). 
In this paper we analyse the implications of the apparent 
disagreement of independent central mass estimates for 
Mkn 501. Starting with an investigation of the constraints 
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given by host galaxy observations (Sect. 2), we proceed 
with an analysis of the robustness of the HE emission es- 
timates (Sect. 3) and finally consider the consequences 
expected from a binary BH modelling (Sect. 4). 

2. Central mass estimates derived from host 
galaxy observations 

Recent mass estimates using host galaxy observations 
mainly rely on the following two statistical correlations: 
(1) Dynamical studies of nearby elliptical galaxies have 
revealed an apparent, almost linear correlation (albeit 
with significant intrinsic scatter) between the central 
BH mass and the B-band luminosity of the bulge 
part of the host galaxy, which is given by Mg ~ 
O.78xlO 8 M (L B , bu i gc /lO lo L 



5) 1 - 08 (Magorrianet al. |1998 



Kormendy & Gebhardt gOOlj , hereafter KG01). 
(2) A much tighter correlation Mjj oc a a seems to exist be- 
tween the BH mass Mh in nearby inactive galaxies and the 
stellar veloc ity dispersion a of the ir hos t bulge (Gebhardt 
Ferrarese & Merritt 200C| ). However, up to 



et al. 2000 



now there is considerable debate over the true slope a. 
Using different samples, Gebhardt et al. (2000) found 
a = 3.75 ± 0.3, while Ferrarese & Merritt ( ^000| ) ob- 
tained a = 4.8 ± 0.54 and a = 4.72 ± 0.36 (Merritt & 



Ferrarese 2001a). Currently, further research is required 
to settle the question whether this difference is mostly 
caused by lower quality data and a less precise regression 
algorithm (cf. Merritt & Ferrarese 2001a) or by systematic 
differences in the velocity dispersions used by the different 



groups for the same galaxies (cf. Tremaine et al. 2002). 
The results using reverberation mapping (RM) (e.g. Kaspi 
et al. |2000t Nelson |200C| ; Wandel |2002[ ) indicate that 
the Mh ~ <J correlation may also hold for nearby AGN. 
However, a critical test of this conclusion depends on both 
a secure measure of the BH mass and an accurate de- 
termination of the stellar velocity dispersion. So far, the 
quality of BH mass estimates from stellar or gas kine- 
matical data (whether ground- or HST-based), which re- 
quire the BH sphere of influence rn — GMh/o 2 to be 
well-resol ved, see ms to increase only modestly (Merritt & 
Ferrarese 2001a ), so that over-estimation may be quite 
possible. Further progress has been expected using RM 
methods (Ferrarese et al. 2001 ). Yet, the accuracy of RM 
may be strongly affected by systematic errors, e.g. due to 
uncertainties in the geometry and kinematics of the BLR 
or due to an unknown angular radiation pattern of the 
line emission, which may result in a systematic error up 
to at least a factor of 3 (cf. Krolik 2001). Moreover, only 
few accurate measurements of a seem to exist for AGN. 



Ferrarese et al. (2001) have recently analysed six AGN 
with well-determined RM BH masses by a careful mea- 
surement of their velocity dispersions and found a general 
consistency with the Mh — a relation for quiescent galax- 
ies. However, only BH masses below ~ 10 8 Mq have been 
included so far, thus leaving out the high mass end of the 
correlation, and in addition, a large scatter is indicated. 
Besides providing a promising tool for the determination 



of BH masses in AGN, the current uncertainties in the 
correlations should be considered, if one tries to assess its 
implication for individual sources such as Mkn 501. 

In the case of Mkn 501, Barth et al. (2002a) have re- 
cently determined the stellar velocity dispersion from the 
calcium triplet lines to be a = (372 ± 18) km/s (cf. also 
Barth et al. 2002b). Applying the Mh — & relations of 
KG01 and Merritt & Ferrarese (2001a), they derived a 
BH mass for Mkn 501 of M H ^ (0.9 - 3.4) x 10 9 M o . 
This mass estimate was supported by the study of Wu 
et al. ( |2002 ), who estimated the velocity dispersions and 
BH masses from the fundamental plane for ellipticals for 
a large AGN sample including 63 BL Lac objects (but 
not Mkn 501). They derived BH masses up to 10 9 M Q , 
but with a potential error up to a factor of two. In par- 
ticular, inspection of the fit in their Fig. 1 indicates a 
possible BH mass for Mkn 501 of - (4 - 7) X 10 8 M for 
M fl (host) = -24.2 mag (Pursimo et al. |2002|) . The general 



challenge of determining a accurately may be illustrated in 
mo re det ail with reference to the recent work by Falomo et 
al. ( 2002 ), who provided a systematical study of the stellar 
velocity dispersion in seven BL Lacs. Using measurements 
in two spectral ranges, they found a velocity dispersion of 
d = (291 ± 13) km/s for Mkn 501, which is significantly 
lower than the one derived by Barth et al. (2002a). Hence, 
if this value is used instead, the BH mass estimated by 
Barth et al. (2002a) is reduced by up to a factor of three, 
i.e. one obtains M H ~ (3.6 - 10.7) x 10 8 M Q . Additional 
support for such a low a- value in Mkn 501 seems to be 
indicated by the original Faber & Jackson relation, which 
yields a ~ 270 km/s (see Fig. 2 in Falomo et al. 2002) for 
Mfj(host) = —24.2 mag. Future research is needed to test 
whether the discrepancy in a is mainly induced by the 
difference in the method deriving a (direct fitting versus 
Fourier quotient routine). 

As noted above, an additional mass estimate for 
Mkn 501 can also be derived from the Mh — ibuigo cor- 
relation. The reported large uncertainties in this re lation 
have recently been examined by McLure & Dunlop ( 2002 ) 
using R-band luminosities, which are less sensitive to ex- 
tinction. By analysing the virial BH masses for a sample 
of 72 AGN, they found the scatter to be quite smaller than 
previously estimated and stressed its usefulness. For ap- 
plication to Mkn 501, we may exploit the absolute R-band 
luminosity of its host galaxy recently derived by Pursimo 
et al. ( 2002 ) (see also Nilsson et al. 1999 ). Assuming 
Hq = 50 km s" 1 Mpc -1 , they obtained Mi?.(host) = -24.2 
mag. If we convert R- to B-band luminosity assuming 
B — R = 1.56 (e.g. Goudfrooij et al 
al I 

L (M B (host)) 



- R = 1.56 (e.g. Goudfrooij et al. |1994| ; Fukugita et 
1995| ; Urry et al. |2000D , we have M B (host) = -22.46 

1.406 x 10 n L r . 



mag, which results in nyiviB^ uosljj = i.wu x iu~i^0. 
Using the KGOl-relation for the B-band luminosity, the 
expected BH mass in Mkn 501 is M H - 1.3 x 10 9 M o , 
but with a potential error of up to at least a factor of 



three. Using the more recent McLure & Dunlop (2002)- 
relation \og(M H /M Q ) = -0.5M fi - 2.96(±0.48), one finds 
M H ^ (0.46 - 4.2) x 1O 9 M . Uncertainties in the de- 
termination of Mfj may further reduce the expected BH 
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mass. For example, values from the literature presented in 
Table 4 of Nilsson et al. ( |1999| ) indicate that M R might be 
up to 0.4 mag higher and therefore Mjj correspondingly 
smaller. More importantly, if a Hubble constant H = 70 
km s _1 Mpc _1 is assumed, one finds Mr = —23.47, which 
results in M H ~ 7.8 x 10 8 M Q (KG01), again with sub- 
stantial scatter of up to at least a factor of three, or 
(2 - 18) x 10 8 M Q (McLure & Dunlop j2002|), thus 



H 



M, 

allowing for a central mass as low as 2 x 10 s Mq 



3. Mass estimates derived in the context of high 
energy emission models 

Following a quite different approach, the BH mass of 
Mkn 501 could also be estimated in the framework of high 
energy emission models: 

(1) With respect to the high energy emission, Fan et 
al. ( 1999 ) have recently determined the central black hole 
masses for several 7-ray loud blazars (including Mkn 501) 
by assuming that the observed 7-rays are produced at 
~ 200 R g . Accordingly, the central BH mass obeys the 
relation Mh/Mq ~ 500 S At, where At denotes the dou- 
bling timescale in seconds and 8 the Doppler factor. For 
Mkn 501 they obtained a central BH mass of Mjj ~ 
0.9 x 1O 7 M , using At = 6 hrs and S ~ 0.9. However, 
current evidence indicates that the relevant timescale At 
might be substantially smaller. Recent observations reveal 
considerable sub-hour variability (e.g. Ghosh et al. 2000; 
Catanese & Sambruna |2000 ) on a timescale as low as 
At = 1200 s (Sambruna et al. 200C). Such a low value for 
the observed timescale could possibly be accommodated 
by assuming a high doppler boosting factor S ~ 18. Hence, 
it appears that the more crucial point in this derivation 
is the assumption that the 7-rays dominating the emis- 
sion are produced at ~ 200 Schwarzschild radii. Indeed, 
at least in the case of the TeV-blazars, the variable, high 
energy emission is usually regarded as produced by mov- 
ing knots or shocks in the j et far from the accretion disk 
(for a review, cf. Mannheim 1997; Aharonian & Volk 2001 



and references therein). While instructive, the derived es- 
timate should thus not be considered as a robust, general 
limit. 

(2) A further mass estim ate for Mkn 501 has been derived 
by DeJager et al. ( 1999 ) following a n app roach developed 
by Hayashida & Miyamoto et al. ( 1998 ). Assuming the 
variation in the accretion process to drive the X-ray and 
TeV variation in the jet via the dynamo effect, their result 
yields a central BH mass of M H = (1 - 6) X 10 7 A/ o for 

5 = (10 — 15). However, due to the absence of a physi- 
cal basis for the required scaling of the Fourier spectrum 
and due to the assumption of a linear proportionality be- 
tween variability timescale and BH mass (cf. Hayashida 

6 Miyamoto et al. 1998), which is probably not valid for 
the blazar class (Kataoka et al. 2001), this estimate again 
does not appear to be robust. 

(3) The estimates (1) and (2) which suggest a BH mass 
less than ~ 6 x 1O 7 M , are strongly model-dependent 
as shown above. We may illustrate this in more detail 



by comparing them with results derived in the context 
of another, high energy emission model. Bednarek, Kirk 
& Mastichiadis (1996) for example, have developed a 
special model for the origin of the high energy parti- 
cles in TeV blazars like Mkn 421 and Mkn 501, assum- 
ing the electrons responsible for the high energy emis- 
sion to be accelerated rectilinearly in an electric field. 
In this model, the mass of the central BH is expected 

to be limited by M H £ 10 8 M Q (£ max /34 TeV) 2 - 5 /4dd' 
where £ max denotes the maximum photon energy and 
^Edd the disk luminosity in units of the Eddington lu- 
minosity. There is strong evidence for a sub-Eddington 
accretion mode in BL Lacs in general (e.g. Cavaliere & 
D'Elia 2002) and particularly for the TeV emitting blazars 
(e.g. Celotti et al. |1998 ). Thus, using characteristic val- 
ues, i.e. Z E dd = ( 0.01 - 0.001) and = 20 TeV (cf. 
Samuelson et al. [1998 ; Konopelko 199E), we arrive at a 
mass M H ~ (2.46 - 14.66) x 10 8 M Q , which is up to ten 
times larger than the estimates (1) and (2). 



4. Mass estimates derived in the context of a 
binary black hole scenario 



In a recent contr ibution (Rieger & Mannheim 2000| , here- 
after RMI; 2001 ), we have shown that the periodicity of 
P bs = 23 days, observed during the 1997 high state of 
Mkn 501, could be plausibly related to the orbital mo- 
tion in a BBHS, provided the jet, which dominates the 
observed emission, emerges from the less massive (sec- 
ondary) BH. If such an interpretation (henceforth called 
the standard scenario) is appropriate, we may derive a 
third estimate for the central mass in Mkn 501. We 
have demonstrated in RMI that, due to relativistic ef- 
fects, the observed period appears drastically shortened, 
so that for the intrinsic keplerian orbital period one finds 
Pk = (6 — 14) yrs . Taking into account that the observed 
emission is periodically modulated by differential doppler 
boosting due to the orbital motion, one may derive a sim- 
ple equation for the required mass dependence in the stan- 
dard scenario (cf. RMI, Eq. (8)): 



M 



D l/3 
obs 



(to + M) 2 /3 (2tt [1 + z] G)!/3 sm « 

yl/(3+a) _ 1 



/1 /(3 + «) + 1 ( 1 -7 cos * ) 2/3 ' 



(1) 



Here a denotes the spectral flux index, z the redshift, v z 
the outflow velocity, i ~ l/T^ the angle of the jet axis to 
the line of sight, ~ 10— 15 (e.g. Mannheim et al. 1996; 
Spada et al. 1999) the bulk Lorentz factor and / the ob- 
served flux ratio between maximum and minimum (for the 
TeV range / ~ 8), while to and M denote the masses of 
the smaller and larger BH, respectively. 
In order to break the degeneracy in this mass ratio, we 
may utilize an additional constraint by assuming that the 
current binary separation d corresponds to the separation 
at which gravitational radiation becomes dominant (cf. 
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RMI). Such a constraint yields an upper limit for the al- 
lowed binary masses and might be associated with the 
key aspect that BL Lac objects are old, more evolved 
and under luminous sources, i.e. they might be close bi- 
naries, probably settled above or near the critical gravita- 
tional separation, because the possibility of removing fur- 
ther angular momentum has been almost terminated as a 
result of declining gas accretion rates. We can specify the 
corresponding gravitational separation d g by equating the 
timcscalc T grav = \d/d\ = (5c 5 /64 G 3 )d 4 /(M m [m + M]) 
on which gravitational radiation shrinks the binary orbit, 
with the dynamical timescale r gas for gas accretion (cf. 



Begelman et al. 1980; note that compared with RMI, this 
estimate for r grav i s a factor 2.5 more precise, cf. Rieger 
& Mannheim 2001). A characteristic measure for r gas is 
given by the Eddington limit r gas ~ 3.77 x 10 7 (77/O.I) 
yrs, assuming a canonical 10% conversion efficiency (cf. 
Krolik 1999p . Using T grav and r gas , one finally arrives at 



3.50 x 10 16 M 8 1/4 m' /4 



m s + Mg) 1 / 4 



(2) 



If we combine this expression with the relevant expression 
for the binary separation (i.e. Eq. (7) of RMI), we obtain 
a second constraint on the allowed mass ratio given by 



(mg + M 8 ) 3 / 4 

,,5/4 1/4 



7.3 x 10' 



/V(3+a) + 1 (1 + z ) 
/VO+a) _ 1 



(3) 



assuming the jet to arise from the less massive BH (for 
the reverse case the BH masses should be interchanged in 
Eq. (Q) and Eq. @). By using Eq. ([[]) and Eq. (g), we may 
determine the binary masses permitted for a BBHS with 
separation near that for which gravitational radiation be- 
comes dominant. In Fig. |l| we have plotted the allowed 
range for the cas e of T & = 15 and a ~ (1.2 — 1.7) (e.g. 
Aharonian et al. |1999|) . The run of the curves is deter- 
mined by Eq. (Q) with the upper bounds given by Eq. ([|). 
We note that the upper bounds depend on the presumed 
timescale for the binary evolution, i.e. on the estimate 
of the gravitational separation. In particular, the derived 
upper bounds may be larger, if, for example, disk-driven 
migration (e.g. Armitage & Natarajan 2002 ) might be rel- 
evant. Fig. [l] reveals that (a) the typical combined mass 
for a BBHS is expected to be less than ~3x 10 8 M©, and 
that (b) a BBHS with, e.g. m~6x 10 7 M© and combined 
central mass (M + m) ~ (1.5 — 2) x 1O 8 M0 appears well- 
conceivable given the current limits. 

On the other hand, if a high central mass of ~ 10 9 M© 
will be established by further research, the proposed bi- 
nary scenario appears to be ruled out. We note, however, 
that even in this case a binary scenario may be still pos- 
sible provided that the jet, which dominates the emission, 
is produced by the primary BH. To illustrate the impli- 
cations in this case, let us consider a (combined) central 
mass of ~ 0.9 x x 10 9 M© by demanding the primary to be 
in the range (5 - 6) x 10 8 M© (see Table 0). The mass of 
the secondary BH then is determined by Eq. (|l|) with the 
masses interchanged, the current separation d by Eq. (3) of 



P ob „=23 days 

r h = 15 



Log(m/1O«M ) 

Fig. 1. Allowed central mass dependence as a function of 
the secondary mass assuming a BBHS, where the observed 
periodicity of 23 days is related to the orbital motion of a 
jet, which emerges from the less massive BH. The allowed 
mass range lies inside the curves and is calculated using 

r b = 15. 

RMI (again with the masses interchanged) and the grav- 
itational separation d g by Eq. The results as shown 
in Table [lj imply a close binary system with d ~ d g . If the 
optically bright QSO stage thus occurs during the binary 
evolution and the applied doppler factors are considered as 
typical, phases of super-Eddington accretion and/or with 
decreased conversion efficiency seem to be necessary for 
the binary to be above its gravitational separation. 



i = 1/Tj, 


1/10 


1/15 


m[10 s M Q ] 


4.22 (3.62) 


3.85 (3.34) 


M [10 s M Q ] 


5.00 (5.00) 


6.00 (6.00) 


d [10 lb cm] 


4.87 (4.76) 


8.54 (8.39) 


d g [10 Ib cm] 


13.1 (12.4) 


13.6 (12.9) 



Table 1. Binary masses for a central core mass of ~ 
0.9 x 10 9 Af©, currently expected separation d and cor- 
responding gravitational separation d g . The calculations 
have been done using a = 1.2 (in brackets: a — 1.7). The 
jet is assumed to be produced by the more massive BH. 



5. Discussion and conclusions 

In this paper, we have considered the apparent discrepan- 
cies between independent central mass estimates for the 
paradigmatic object Mkn 501. Our main conclusions can 
be summarized as follows: 

(1) Estimates inferred from host galaxy observations using 
the reported, statistical correlations for the central veloc- 
ity dispersion and the bulge luminosity reveal a tendency 
for the central BH mass in Mkn 501 to be larger than 
~ 5 x 10 8 Af©. If verified by further investigations, this 
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would exclude both the standard binary scenario, where 
the jet is emitted from the less massive BH, and several 
high energy emission models proposed so far. However, as 
the correlations still show substantial intrinsic scatter and 
large uncertainties due to the use of limited samples - es- 
pecially with respect to the high mass end - the accuracy 
of this estimate should be treated with caution. In partic- 
ular, we have shown that the central mass limit could be 
easily up to a factor of (2 — 3) smaller, thus allowing for 
the standard binary scenario, in which case a simple ex- 
planation for the divergent mass estimates from emission 
models and host observations appears possible. 

(2) Independent mass estimates from high energy emission 
models suggest a mass for the jet emitting BH in Mkn 501 
which is smaller than ~ 6 x 10 7 Af Q . However, this esti- 
mate is quite model-dependent and so cannot be used as a 
universally valid constraint. In particular, emission mod- 
els have been developed for Mkn 501 where a central BH 
mass up to ten times larger appears quite possible. 

(3) If a binary BH interpretation for the periodical vari- 
ability (cf. RMI) is appropriate in the case of Mkn 501, 
the jet dominating the emission has to be produced by the 
less massive (secondary) BH. Using characteristic jet pa- 
rameters and assuming both Eddington-limited accretion 
and a binary separation comparable to the gravitational 
one, we have shown that the combined central (primary + 
secondary BH) mass should be smaller than ~ 3 x 10 8 M©. 
The binary model particularly permits a system with 
m = 6 x 10 7 Af© and (m + M) ~ (1.5 - 2) x 10 8 M Q , 
possibly allowing for a convergence of constraints from 
emission models and host observations. 

(4) If further research clearly establishes a high central 
BH mass ~ 10 9 M© in Mkn 501, this would call for a 
change of fundamental assumptions and/or a modifica- 
tion of the parameter space regarded to be typical in the 
context of several high energy emission models. The only 
way to avoid this conclusion seems to be the assumption of 
a highly unequal BBHS with the jet dominating the emis- 
sion produced by the less massive and the central mass 
dominated by the primary BH. However, it then appears 
no longer be possible to explain the observed periodicity 
via the orbital motion. On the other hand, even in the case 
of a high central mass, a binary model for the observed 
periodicity may be still possible provided the jet is pro- 
duced by the primary BH. Yet, for the binary to be near 
or above the separation at which gravitational radiation 
becomes dominant, one then requires the earlier binary 
evolution to pass through phases of super-Eddington ac- 
cretion and/or decreased conversion efficiency r\ < 0.1. 
Such conditions need not necessarily be ad hoc but have 
already been considered in the context of galaxy evolution 
and the growth of massive BH (e.g. Collin et al. 2002; Yu 
& Tremaine |2002] ). 

In view of Mkn 501 as a paradigmatic object, bringing 
its mass estimates to convergence continues to be an im- 
portant task. Further research on the parameter space of 
emission models and the careful observations of nearby 
galaxies will be particularly valuable in order to minimise 



the uncertainties in the statistical correlations, to evaluate 
the impact of host observations on high energy emission 
models, and to assess the plausibility of a BBH model. 
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